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Abstract 

We study the late time evolution of flat and negatively curved Friedmann- 
Robertson- Walker (FRW) models with a perfect fluid matter source and a scalar 
field arising in the conformal frame of f{R) theories nonminimally coupled to 
matter. Under mild assumptions on the potential V we prove that equilibria cor- 
responding to non-negative local minima for V are asymptotically stable, as well 
as horizontal asymptotes approached from above by V. We classify all cases of 
the flat model where one of the matter components eventually dominates. In par- 
ticular for a nondegenerate minimum of the potential with zero critical value we 
prove in detail that if 7, the parameter of the equation of state is larger than one, 
then there is a transfer of energy from the fluid to the scalar field and the later 
eventually dominates in a generic way. 



1 Introduction 



In cosmological models containing scalar fields the exponential potential function, V (0) = 
Voe"''*'^, is the most popular not only because of the variety of alternative theories of grav- 
ity which predict exponential potentials, but also due to the fact that this potential has 
the nice property that V (x V which allows for the introduction of normalized variables 
according to the formalism of Wainwright et al [Q. Among many investigations with 
exponential potentials, we mention [|13] and references therein for FRW models; for 



spatially homogeneous Bianchi cosmologies with an exponential potential; and |jT7 
for models containing both a perfect fluid of ordinary matter and a scalar field with 
an exponential potential, the so-called "scaling" cosmologies and for nonminimally 
coupled scalar field; and references therein for scalar-tensor theories with exponen- 
tial potential; ||22| for inhomogeneous cosmologies; |jl5| for flat and open FRW models 
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and [0] for Bianchi type I models with multiple fields containing as a special case that 
of a multi-exponential potential studied in the context of assisted inflation; |jl6| for an 
elegant mathematical generalization of multi-exponential potentials; |^ for an inves- 
tigation of a class of exponential potentials, allowing exact solutions in the context of 
dark energy; for scalar-tensor phantom cosmologies with exponential potentials; H] 



for the reconstruction of the f{R) theory starting from the exponential potential in an 
effort to solve the puzzle of cosmic acceleration. However, a large class of potentials 
used in scalar-field cosmological models have a local minimum. Examples of potentials 
belonging to this class are polynomial potentials of the form V (0) = A0^", generalised 
and logarithmic potentials V (0) oc 0" (In^)™' or V (0) oc 0"'exp (— A0™') studied in 
potentials used in the study of isotropisation of Bianchi type models in scalar-tensor 

potentials V (0) oc Ae~'^'^ + Be~^'^ in an effort to avoid eternal ac- 
in scalar-tensor quintessence (see also for potentials reconstructed 



theories and pT] 



celeration 48 



from observations), in phantom cosmology ||T0[, in double scalar-tensor cosmologies [0. 
Other important examples include chameleon effective potentials and potentials in 
conformally related theories of gravity, for example 



2/3<}!.y 



(1) 



which arises in the conformal frame of the R+aB? theory ||3^, Multivalued potentials 
arising in realistic f{R) gravity theories were studied in 
potentials may drive a flat FRW universe to recollapse [p7 



Unbounded from below 



Since the nature of the scalar field supposed to cause accelerated expansion is un- 
known, it is important to investigate the general properties shared by all FRW models 
with a scalar field irrespective of the particular choice of the potential. Exact solutions 
for fiat FRW models containing only a scalar field with arbitrary potential were obtained 
in (see also for Bianchi type I and V models containing ordinary matter). Nev- 
ertheless, the number of papers with mathematically rigorous results is small. Flat FRW 
models having an arbitrary potential with a positive lower bound were studied by Foster 
Generalisation to models including ordinary matter and scalar fields with potentials 

3|1 studied the asymptotic 
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having a zero local minimum was presented in Rendall 
behaviour of homogeneous models with a scalar field having an arbitrary potential with 
a positive lower bound and showed that the no-hair theorem holds for that case (see also 
Bianchi type I- VIII models were shown to isotropise in the context of k-essence 
51. The oscillatory behaviour of scalar fields with a general potential of the form 
0^/2 + (0^) and ordinary matter was investigated in |^ for fiat FRW models. 

and the 
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m 

V(0) 
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Collapsing models were built using homogeneous scalar field solutions in 
case where a scalar field is coupled to a perfect fluid was studied in . 

The motivation of this investigation comes from a number of physical theories which 
predict the presence of a scalar field coupled to matter. For example, in string theory 
the dilaton field is generally coupled to matter Nonminimally coupling occurs 

also in scalar-tensor theories of gravity in higher order gravity (HOG) theories 
and in models of chameleon gravity [^. In particular, for HOG theories derived from 
Lagrangians of the form 

L = f{R)^ + 2L^m, (2) 
it is well known that under the conformal transformation, 'g^y = f (R) g^^, the field 
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equations reduce to the Einstein field equations with a scalar field (j) as an additional 
matter source, where 

l^^fiR)- (3) 

We do not enter into the discussion about the regularity of the conformal transformation, 
or the equivalence issue of the two frames (see for example |jl9[ and references 



therein). Assuming that (|]) can be solved for R to obtain a function R (0) , the potential 
of the scalar field is given by 

ViR{<P)) = ^{Rf-f), (4) 

and quadratic gravity with the potential (|ip is a typical example. The restrictions on 
the potential in the papers [|^, 0, ^ were used in |^ to impose conditions on 
the function f {R) with corresponding potential The conformal equivalence can 
be formally obtained by conformally transforming the Lagrangian (|^) and the resulting 
action becomes IQ, 

S = J d^x^ ^R - [{d<pf + 2V (0)] + 26^2^273*/:^ [e-^^'^g, ^) } . (5) 

Variation of S with respect to 'g yields the field equations, 

G^. = T^.(5^,0) + T^.(5^,^), (6) 
and variation of S with respect to (j) yields the equation of motion of the scalar field. 

Note that the Bianchi identities imply that 

V^T;,(fi^,vI/)^0, V^T^,(5^,0)^O, (8) 

and therefore there is an energy exchange between the scalar field and ordinary matter. 

In this paper we study the late time evolution of initially expanding fiat and negatively 
curved FRW models with a scalar field having an arbitrary bounded from below potential 
function V (0) . Ordinary matter is described by a barotropic fiuid with equation of state 

p=(7-l)p, 0<7<2. 

The scalar field is nonminimally coupled to matter according to (j^), (^) and (^. Under 
general assumptions on the potential function V{(f)) (see Assumption ffl) we first show 



in Section |2.1| that stable solutions are related to (possibly degenerate) local minima 
of V with non-negative critical value and in cases where V approaches a horizontal 
asymptote from above. In Section ^ we focus for the sake of simplicity on the fiat case 
and investigate the energy exchange between the fiuid and the scalar field in the case of 
the exponential potential and in the case of a nondegenerate local minimum of V . In 
the former case (Section |3.1| ) the fiuid energy never dominates, nevertheless the scalar 
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energy always contributes as a nontrivial fraction of the total energy and in some cases 
it totally dominates over the fluid, but stable equilibria may also appear corresponding 
to scaling solutions, similarly to the uncoupled case treated in Concerning the case 
of nondegenerate minimum of V (Section |3.2|) we perform a careful study of the system 
and prove in detail that the late time qualitative behaviour is determined by the sign 
of the pressure of the fluid. In particular, the energy exchange between the two matter 
components is such that for 7 > 1 the scalar field eventually dominates (except possibly 
for a particular solution, see Theorem 0). To our knowledge there is no rigorous proof 
in the literature of this result and even for the uncoupled case studied in Section 4] 
the same conclusion that we find here, though reasonable, was only conjectured and no 
proof was provided. 



2 Flat and negatively curved FRW with an arbitrary 
potential 

For homogeneous and isotropic spacetimes the field equations (^ reduce to the Fried- 
mann equation. 



and the Raychaudhuri equation, 

H^-li^-l,.^,, (10) 
while the equation of motion of the scalar field (|^), becomes 

;p + sH<p + v'{<p) = ^^p. (11) 

The Bianchi identities yield the conservation equation, 

p + 3jpH = -^^p<P, (12) 
(see for example Q). For simplicity we drop the tilde from all quantities. We adopt 



the metric and curvature conventions of [50|. a (t) is the scale factor, an overdot denotes 



differentiation with respect to time t, H = a/ a and units have been chosen so that 
c = 1 = SttG. Here V (0) is the potential energy of the scalar field and V = dV/dcj). 

In the following, we are going to study the late time behaviour of solutions of (PH])- 
(p!^), which are expanding at some initial time of observation, i.e. H{Q) > 0. For fiat, 

k = 0, models the state vector of the system (|T0|)-(|T2|) is ^0, 0, p, ifj, i.e. we have a 

four-dimensional dynamical system subject to the constraint (^. Defining y := (p and 
setting (4 — 37)/v^ =: a, we write the autonomous system as 

= 1/, 

y = -3Hy - V (0) + ap, 

p = -3-fpH - apy, (13) 

TT 1 2 7 
^ = "2^ -2^' 
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subject to the constraint 

m' = p+ly^ + V{<P). (14) 

Remark 1 The function W defined by 

Ty (0, y, p,H) = H'-^- + y^^^ + ^ (15) 

satisfies 

W = -2HW. (16) 

By standard arguments in ordinary differential equations theory it follows that sgn{W) 
is invariant under the flow of (|T3p . We deduce that solutions with W positive, null, or 
negative, represent scalar field cosmologies with k = —1,0, 1 respectively. 

Remark 2 Similar arguments applied to the third of ( |I3D show that if p > at some 
initial time to, then p{t) > throughout the solution. Notice that in scalar field cos- 
mologies in the context of GR, the right hand sides (rhs) of (IJ) and ^TSj ) are zero. The 
energy exchange between the two matter components (cf is reflected to the time 
derivative of the energy density of the scalar field, 

which contains the extra term apy. Furthermore, 

e + p=-3H{y' + jp), (17) 
implies that, for expanding models the total energy e + p of the system decreases. 



2.1 Stable equilibria for open spatial topologies 

The equihbria of (p!3|) are given by (0 = 0*,?/ = 0,p = 0,H = ±a/V^(0*)/3) where 
V'{(f)^:) = 0, and we are going to discuss their stabihty for expanding cosmologies {H > 0) 
with open spatial topology {k = — 1,O).0 Hereafter we are dealing with a quite general 
class of potentials which includes, for instance, polynomial functions with even leading 
term and exponential-decay functions. Essentially all we need is that V{(j)) is eventually 
non-negative as — )■ ±oo and has a finite number of critical points: 

Assumption 1 We suppose that V{(j)) G C^(]R) satisfies the following conditions: 

1. the (possibly empty) set {0 : V^(0) < 0} is bounded; 

2. the (possibly empty) set of critical points ofV{(j)) is finite. 

^Thc function V is said to have a degenerate local minimum at 0, if both V' and V" vanish at 0, . 
Moreover, if y, p — >■ as t +oo and if) — > +oo, while limt_).+oo -ff — > ± lim0_i._|_oo y/V{(j))/3), then we 
call the asymptotic state {(f> — +oo,y = 0, p = 0, H = ±-\/y(cx))/3), stable equilibrium "at infinity", 
although a more appropriate term should be stable configuration. 
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Critical points of V with negative critical value are not equilibria and they rather 
allow for recoUapse of the model. Moreover, nondegenerate maximum points (with non 
negative critical value) for V are unstable, as can be easily seen by linearizing system 
(0) at the corresponding equilibria and verifying the existence of at least one eigenvalue 
with positive real part. If the maximum points for V are degenerate, it can be shown that 
they are also unstable. In fact, the function e + p — is indefinite in the intersection 

of the set {VT > 0} with an arbitrary neighborhood of the equilibrium point and strictly 
decreasing under the flow of the system for expanding cosmologies (see eq. (|l^). On the 
contrary, it can be proved that equilibria corresponding to local minima of the potential 
with non-negative critical value are asymptotically stable. More precisely, we have the 
following result. 

Proposition 1 Let (j)^ a strict local minimum for V{(j)), possibly degenerate, with non- 
negative critical value. Then, y** = (0*, 1/* = 0, p* = 0, = \J~^^^) is an asymptotically 
stable equilibrium point for expanding cosmologies in the open spatial topologies k = 
and k = —1. 



The proof of this Proposition is given in the Appendix. 

Of course, we must also consider possible stable configurations corresponding to a 
diverging value for the field 0. In view of the Assumptions, this can happen only when 
V{(j)) has critical points "at infinity" with asymptotic value i > 0. For simplicity let 
us restrict to the case (p ^ +oo. Then it can be shown that when V{(j)) — )■ i~ , the 
configurations is unstable. The proof can be done by introducing a new scalar coordinate 
ip = V{(j)) — £, and writing the system in terms of ip. In this way the critical point at 
infinity becomes a finite local maximum {ip = 0). Arguing similarly as sketched before 
for local maxima, it is shown that when V{(p) — )■ i~ , the configuration is unstable. 
On the other side, when V{(p) — )■ (e.g. as it happens for the exponential potential 
V {(p) = Voe~^^) the critical point at infinity is asymptotically stable, i.e., ?/,p — )■ 0, and 
(p — i- +00 as t — > +00. Moreover, since lim(^_^+oo V (0) = we conclude from ( p!4D that 
also H ^/i+JS as t — cxD. The proof follows after suitable adaptation of the arguments 
used in Proposition 4 in |ET|. 



3 Transition from the matter to the scalar phase 

In the following we study the energy transfer from the perfect fluid to the scalar field. 
We are interested to study the late time behaviour near asymptotically stable equilib- 
rium configurations. We suppose that the initial data in the basin of attraction of this 
equilibrium are such that the fluid is the dominant matter component, i.e. 

Po > eo, 

and we are asking whether there is a time ti such that 

e{t)>p{t), Vt>ti. (18) 

This question is closely related to the classification of scalar field potentials with cosmo- 
logical scaling solutions [^. It is also relevant to the coincidence problem, that is, why 
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dark energy and matter appear to have roughly the same energy density today (see 
and references therein). 

There are two cases to be examined for the hmiting value of V{(j)), regardless of the 
critical point being finite or "at infinity" : 

(i) V{(f)) asymptotically tends to a strictly positive limit. Then, it is easily seen that 
the transition (|18D does always happen, since the critical value of the potential behaves 
as an effective cosmological constant and the energy of the scalar field tends to this value 
whereas the energy of the fluid tends to zero.| 

(ii) V{(f)) asymptotically tends to zero. This is the nontrivial case and from now 
on we will focus ourselves to this question, studying both the exponential potential 
and nondegenerate minima. For simplicity we will restrict ourselves to the fiat, k = 0, 
case, although it seems that the same results can be stated for the general open spatial 
topology case []3T|. 



3.1 Exponential potential 

In the case of the exponential potential V{(f)) = Voe~^'^, we already know that — )■ +oo 



and ?/,p — )■ 0, (see the remarks at the end of Section |2Tl|). To investigate the late time 



behaviour of the energy we follow the ideas from [17], writing the system in the variables 



'K)e-t<A y 

w ' ' 




3 H ' y^H 

and H, with a new time variable r defined by dr = 3Hdt: 

-— = w [ := + z H — (1 — w — z 

dr V 76 2^ 



^ = z(-l + z' + ^{l-w'- z')) + ^ (Xw' + a{l -w'- z')) , (19) 
dr V 2 / v6 

^ rr(.2 . 7 

dr " 



h[z^ + 1{1-w^-z^)), 



where we have used the constraint (0) to eliminate p. This system has the nice prop- 
erty that the third equation decouples from the other two. Obviously the scalar field 
dominates when w"^ + z"^ — )■ 1. Studying the Jacobian at the equilibria of (p!9[) one can 
find that there are three possible stable configurations, depending on the couple (A, 7). 
A summary picture is sketched in Figure |l]. In particular: 

1. {w, z) = 1 — is a scalar field dominated solution and is the only stable 

asymptotic configuration when 7 < ^^i^^! 



^In a flat FRW model where the perfect fluid at early times is that of radiation and later that of 
matter, dark energy domination during the radiation epoch should be avoided. Therefore, the relevant 
question is when the scalar dominance would start. Unfortunately, qualitative methods as those used 
in this study are unable to answer this kind of questions. 
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Figure 1: Summary of the different stable configurations when V{(f)) = Vqc The 
three regions depicted here corresponds to the three cases hsted on page ^ 



(w, z) 



^0, ^(2^-^) ) ^ scahng solutions, where neither the scalar nor the fluid 
energy eventually dominates, and is the only stable asymptotic configuration when 

X ^ 9 /2 ^ 4(4-^/6A) . 

3. {w,z) = ' , _4_)_3^Tt-^A ) another scaling solution that is the 

only stable asymptotic configuration when (A, 7) is on the complementary of the 
other two regions. 

We can therefore conclude that, in the case of the equilibrium "at infinity" determined 
by an exponential-type potential, the scalar energy e generically remains as a nonzero 
fraction of the total energy. In the first case e totally dominates over the fluid energy p, 
in the other cases both energies give a nontrivial contribution to the total energy (see also 



32| for scaling solutions). Interestingly enough, these are the same conclusions of in 
the case of a scalar field non coupled to matter, although the ranges of the parameters 
(A, 7) are different. 



3.2 Nondegenerate local minima of V{(j)) 

Without loss of generality we suppose that 0* = and therefore, the potential can be 
written in a neighborhood of = as 

V^(0) = ^AV + W), A>0. (20) 

It must be stressed that this is the most general form of an arbitrary potential near the 
equilibrium 0=,,. In fact the first assumption does not enter in the study of the late time 
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behaviour around a critical point 0* G M, because for that situation only the behaviour 
of the potential near the point is important and no further assumptions on V are actually 
needed. For the sake of simplicity we will neglect higher order terms in the potential 
V{(l)). Nevertheless, if we take them into account, it can be shown by lengthy calculations 
that the results we are going to state are not affected. As done in the previous case we 
introduce new variables w, z 



A I y .^.s 



and rewrite the system as 



w = Xz + 3Hw ( + ?(1 -w^ ~z^ 



-\w + 3Hz{z^ + ^{l-w^- + 3H \^-z + ^{1 - w'' - z^) ] , (22) 



H = -3H' ^z^ + l(l-w^-z^' 

which unfortunately does not possess the decoupling property of the corresponding sys- 
tem (|I9|). Furthermore, we cannot use the time variable r unless we lose regularity prop- 
erties of the system and we cannot infer about the stability of the equilibria by looking 



at the Jacobian of (p2D. Actually we will tackle this nontrivial situation performing a 



qualitative study of (]22|) and prove the following result. 

Theorem 2 Let 0^, be a nondegenerate minimum ofV{(f)) with zero critical value. Con- 
sider the solutions of (|1^) with k = approaching the (asymptotically stable) equilibrium 
point (0=,,, y = 0, p = 0, H = 0) . Then if ■y < 1, for every such solution the fluid energy p 
eventually dominates over the scalar energy e, whereas z/7 > 1, e eventually dominates 
over p in a generic way, i.e., except at most for a particular solution of the system. 

Proof. As said before, suppose without loss of generality that V^cj)) is as in (^), and 
consider the variables (|2T|) that bring the system in the form (^21). We first define the 
(positive) function Rhy = w"^ + z^ (of course initial data such that -R(O) < 1 imphes 
Rit) < 1, Vt > 0) and divide the proof into some steps. 

(1) If '-f < 1 then i? — 7- for every solution. 

Let us observe that the flow of the system eventually lives in the compact set 
{{w,z,H) : w"^ + z^ < 1, H ^ [0,Hq]} and that the w-limit point of every solution, 
is such that H = 0. Therefore, if (wq, Zq, 0) is a limit point then it is easily seen by ( p2D 



that the circle w"^ + z"^ = Wq + Zq, H = is made by limit points, which means that the 
solution asymptotically approaches this circle. Suppose by contradiction the existence of 
a solution such that Wq + Zq> 0. Then the only possibility is that Wq + Zq = 1. Indeed, 
passing to polar coordinates 

w = R cos 9, z = RsinO 
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we obtain 



e + X = 3Hcose { -sme+ ^ ^ ^ 



R = 3H(l - R^] 



V6 R 

RC^-sin'^e) +-^sm0 
V2 J 



(23) 



H = -3H^ i R^ sin^ O+'^il-R^' 



If the solution approaches R = Rq > then the first equation says that 6 asymptotically 
behaves like —At. This fact implies the following asymptotic estimates as t — )■ +oo (recall 
that H monotonically) : 

^ 3if - sin^ dr = log a^^^'^^/^ , 3H sin 6 dr < +oo (24) 

and therefore dividing both sides of the second equation in (|23| ) by (1 — R^) and using 
R ^ Rq, the second equation can be integrated to give 

hm a(t)3(^-i)«o - 



t-S> + 00 1 — Rq 

which is consistent with Rq > only when Rq = 1. But this fact implies 7 > 1 (when 
7 = 1 the argument cannot be applied, so we must include this possibility, see Remark 
^ below), a contradiction. 

(2) If '-f > 1 and i? — )■ then R/H is eventually hounded. 

Let us begin assuming H/R — )■ 0. Then, as above, 9 = —At and estimates ( p4|) hold. 
Using — )■ in the third equation of ( p3D we see that H{t) = 2/3'jt and therefore 



Using again _R — )■ we can neglect superlinear contribution for R in (^3]) , and plugging 
in the above estimate we obtain 

R(t)^(^^^ ^ R{to) + c s^~~ sin e{s)ds 

for an appropriate constant ci. Possibly shifting the initial instant to we can assume that 
the quantity in square brackets above is strictly positive, which implies that R(t) — )■ +00 
if 7 > 1, a contradiction. Then H/R -f^ 0. To show that this implies that HjR is 
bounded away from zero, i.e., RjH bounded, we set S = H/R and consider system (p2D 
in the new variables i?, if, 5", X and F, where 

w = ^RX, z = ^RY. 
3a 3a 
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By hypothesis, the cj-hmit points of the solution are such that R = 0, H = and since 
the cj-limit set is invariant, the triple {S, X, Y) of a limit point must belong to a solution 
of the system 

S = -SV, 

X = -Y{-X + SX), 

Y = X{-X + SX), 

as it can be easily seen writing the system in the new variables and setting R = H = 
0. If there exists an w-limit point of the kind {R = 0,H = 0,S = 0,Xo,lo) then 
using the above system we see that every point such that R = 0, H = 0, S = 0, and 
_j_ y2 _ -^2 _|_ is an oj-limit point, that is, the solution curve approaches the circle 
{R = 0,H = 0,S = 0,X, Y) with X^ + Y^ = X^ + Y^, and in particular S ^ 0. 

(3) If •y > 1 there exists at most one solution such that i? — )■ while every other 
solution is such that i? — )■ 1. 

Assume by contradiction the existence of two distinct solutions ^, rj such that R{t) — )■ 
0. By the aforesaid H/R is bounded away from zero (equivalently, R/H is bounded) for 
both ^, f] and therefore we cannot exploit anymore the above argument to show 6 = —At. 
A change of coordinates is then needed in order to control the behaviour of 6 for these 
solutions. In the following we will derive the appropriate change of variables for a general 
system of the form 

ti) = Xz + Fi{w,z,H), (25a) 
z = -Xw + F2{w,z,H), (25b) 
H = F3iw,z,H), (25c) 

where Fi{w,z,H) are C°° functions such that 

Fi(0,0,0) = 0, (Fi),»(0,0,0) = 0, (F,),,(0,0,0) = 0, (26) 



and F3 is a negative function for H > 0. It is easily seen that (p^) has the above form. 
Passing to polar coordinates w = RcosO, z = RsinO we obtain 

■ , F2(Rcose,Rsme,H)cos9 - Fi(Rcose,Rsme,H)sme , , 

e + A = ^ , (27a) 

R = Fi{R cos 9, R sin 9, H) cos 9 + F2{R cos 9, R sin 9, H) sin 9, (27b) 

and (pGl) together with ( p7aD suggests that the condition 

F,(0,0,if) = 0, 2 = 1,2. (28) 

would ensure 9 = —At. Since (^) does not hold a-priori, and in particular does not 
hold for the system we are considering, we look for a coordinate change 

Wi=w-fiH), Zi = z-g{H), 
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where /, g are opportune functions of H. Writing the system in the new unknown 
functions wi,zi,H and imposing (|28D for the new system we obtain 



Xg + F,{f, g, H) - f{H)F^{f, g, H) = 0, (29a) 
-A/ + F2(/, g, H) - g'{H)F,{f, g, H) = 0. (29b) 



This amounts to look for a solution of the original system ( ^5a| )-( p5q) . Indeed, if 
{w{t), z{t), H(t)) is such a solution and denoting by t{H) the inverse function of H{t), 
we see that 

fiH) = witiH)),giH) = zitiH)), 

satisfy (|29aD -( pE|) . 

Coming back to the original problem, we will use the solution rj to perform the 
coordinate change described above and apply it to the other solution ^. Notice also that, 
since i? — )■ and R/H is bounded, it is easy to verify that the conditions 

lim f(H) = lim g(H) = 0, (30a) 
Hf'{H), Hg'{H) bounded near H = 0, (30b) 

hold for f,g. Passing to polar coordinates, the new 6 asymptotically behaves like —Xt, 
and it can be found with lengthy but straightforward calculations that the new R := 
wf + zf satisfies the equation 

R = SHR^^R^ sin ^2 ^ 2^i?sin^ + g^ + '^{1 - R^ - f - g^ - 2/i?cos^ - 2gRsme) + 
[(/ + Hf\H)) cos e+{g + Hg\H)) sin 6] RsinO + 2g sin - |(i? + 2/ cos e + 2g sin 6) 

- sin^ ^ - ^ sin e{R + 2/ cos 6 + 2g sin 6) I . (31) 

Therefore under the assumption that — )■ (note that the new R goes to zero if and 
only if the old R does, since f,g^O) we obtain (using also ( P0a| )- (|30b|) ) 

R = 3HR Q - sin^ 6 

thus, R = a^*-^"^-'/^ which is consistent only if 7 < 1 (again 7 = 1 does not allow to 
perform the above argument and so must be included, see again Remark |^ below), a 
contradiction. Therefore, except r] every other solution is such that R does not tend 
to zero, and a similar argument as in the step (1) ensures that actually — )■ 1. This 
concludes the proof. ■ 



+ 



Remark 3 We note that the argument of the proof relies on estimates of the type 
which are consistent when 7 7^ 1 . Therefore we stress that the above argument fails to 
give information in the case 7 = 1 when actually both situations, t <^ p and p ^ e, may 
generically occur depending on the initial data of the problem. 
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4 Discussion 



We analysed the late time evolution of flat and negatively curved expanding FRW models 
having a scalar field coupled to matter. We proved that, equilibria corresponding to non- 
negative local minima of V are asymptotically stable. In case the minimum is positive, 
say ^(0*) > 0, the energy density, e, of the scalar field eventually rules over the energy 
density of the fluid, p, and the asymptotic state has an effective cosmological constant 
V{(j)^:). In case the minimum is zero and nondegenerate, then p eventually dominates 
over e if 7 < 1 and e dominates over p if 7 > 1. This result could be interesting in 
investigations of cosmological scenarios in which the energy density of the scalar field 
mimics the background energy density. For viable dark energy models, it is necessary 
that the energy density of the scalar field remains insignificant during most of the history 
of the universe and emerges only at late times to account for the current acceleration of 
the universe. 

It must be also stressed that the energy transfer from the scalar field to the fluid 
showed in Theorem ^ is proved generically, i.e. even in the 7 > 1 case there is possibly a 
special solution such that p dominates over e, nevertheless this solution is unstable under 
perturbation of initial data of the problem. The generic behaviour of solutions is not a 
new feature of models involving scalar fields. We mention for example the well known 



instability (with respect to bounded variation initial data) result by Christodoulou |T2 



about scalar fields collapsing to a naked singularity and even in the realm of homoge- 
neous stars (regularly matched with nonhomogeneous exteriors), genericity of collapsing 
solutions ending in a black hole has been discovered to be true both in GR p8| and in 
HOG gravity |g. 



The above results were rigorously proved only assuming that critical points are finite, 
and that V{(f)) is eventually non-negative as — )■ ±00. It must be remarked that the 
latter assumption does not enter in the study of the late time behaviour around a critical 
point 0* G M, because for that situation only the behaviour of the potential near 0* is 
important and no growth at infinity assumptions on V are actually needed. The other 
assumption plays the role to ensure that the flat plateau at infinity are unstable (as 
proved for the special case of the potential (|I]) in |^), whereas if V^cj)) — )■ i from above 
we have an asymptotically stable configuration. As a representative of the latter class of 
stable configurations we have studied the late time behaviour of the energy in the case 



of exponential potential, finding results similar to the non coupled case [|T7|- However 



our dynamical system is much harder to deal with, because the evolution equation for 
H does not decouple as in and moreover its equilibria are not hyperbolic. When 
the decay parameter A is small, where "small" depends on 7, see region 1 in Figure 
|I], the scalar field dominates, whereas other stable configurations correspond to scaling 
solutions where both energies nontrivially contribute to the total energy of the system. 

One can encounter in the same potential all the situations studied so far, for instance, 
one can have a potential exhibiting more than one local minima, or a local minimum and 
an exponential decay at infinity. In that case since at least two asymptotically stable 
configurations for the system can be reached, the late time behaviour of the energies and 
their exchange necessarily depends on the initial data of the problem. 

The cases studied in Section ^ of course do not cover all possible situations, even 
under Assumption |l]. In particular, one can take into account degenerate minima for the 
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potential. In this case, it is not hard to see that the counterpart of systems ( ^2]) and (^31) 
becomes irregular near w = z = 0, and to infer about the behaviour of the angle 6, one 
cannot straightforwardly extend the techniques used here. Nevertheless, research in this 
direction tends to confirm the same generic behaviour showed here for the nondegenerate 
case [^. Therefore, it seems appropriate to conclude that energy transfer is completely 
driven exclusively by the sign of the fluid pressure when the scalar reaches a finite stable 
equilibrium, whereas for stable configurations "at infinity" the potential profile becomes 
crucial. 

Another important question that should be further investigated is the case of closed 
cosmologies. We believe that a closed model cannot avoid recoUapse, unless the minimum 
of the potential is strictly positive. In that case, the asymptotic state must be de Sitter 
space. 
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Appendix 

Proof of Proposition |l]. 

We will follow the ideas used in in the absent of ordinary matter. We start the 
proof by considering the case V^(0*) > 0. Let V > V{(f)^) be a regular value for V such 
that the connected component of — oo,\^]) that contains 0=,,, is a compact set in 

M. Let us denote by A C M this compact set and define fl as 

n = {{4>,y,p,H) : 4>eA, {e + p)<V,p>0, W{^,y,p,H) G [0,W]}, 

where is a positive constant. Then is a compact set because the following facts 
hold: 

1. is a closed set in M^; 

2. V{4>,) < V{(f)) <V,W4>e A; 

3. + V{(j)^) < 1?/^ + V{(f)) + p <V, and therefore y is bounded; 

4. p <V — — V{(j)) <V — V{(f)^:), and therefore p also is bounded; 

5. from ( [T3| ) and the above facts, < <W + V, and so H is bounded too. 
Let f2+ C i7 be the connected component of Q containing p*. It is easy to see that H 



is positive on fi+ and we claim that fi+ is positively invariant with respect to (0), i.e. 
solutions, with initial data in live in Vt > 0. Indeed, let x(t) be such a solution 
and i = sup{t > : H{t) > 0} G MU {+cx)}. When t < t, equations (^) and (|T7p imply 
that both W and e + p decrease. Moreover, it can be proved by contradiction that 

G A, Vt < t, (32) 
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otherwise there would exist some t <t such that V{(f){t) > V, but then 

V < V{<Pit)) < ^y\t) + VicPit) + p{t) = e{t) + p{t) < V, 

a contradiction. Thus, (|32D holds. But since W > along the flow (cf Remark |l]), it 
follows that 



H{ty > -{y{ty + Vim) + pit)) > > 



where the last inequality follows from (|32|). All in all we have proved that, as long as H 
remains positive it is strictly bounded away from zero, thus, t = +oo, and from this fact 
it is straightforward to conclude that x(t) lives in f2+ Vt > 0. 



Using the properties of i7+ proved so far, LaSalle's invariance theorem |£2[ can be 
applied to the functions W and (p + e) in Q^, to show that every solution with initial 
data in f2+ must be such that HW — )■ and Hiy"^ + •yp) — )• as t — !■ +oo. Since H is 
strictly bounded away from zero in Q^, both W and (y^ + 7p) must go to zero, which 
means y -)■ 0, p -)■ and H'^ - lV{(f)) 0. Now the fourth of (|13D 

H = -liy' + 7P)-W, (33) 

implies that H is monotone and therefore admits a limit. This means that V{(j)) also 
admits limit, and this limit must be unavoidably ^(0^,), otherwise V'{(f)) would tend 
to a nonzero value, and so would the right hand side in the second equation of (0), a 
contradiction. Therefore the solution approaches the equilibrium point p*. 

If V{(t>^) = 0, the above argument can be easily adapted. In this case the set Q is 
connected and we choose fi+ to be its subset characterized by the property H >0. The 
only point in with if = is exactly the equilibrium point p*, and so if H{t) — )■ the 
solution is forced to approach the equilibrium since H is monotone; if by contradiction 
H(t) had a strictly positive limit, we could argue as before to flnd ?/ — )■ 0, p — )• and 
ly — 7- and so H must necessarily converge to zero. ■ 
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